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Abstract 

How to give a natural geometric definition of a covariant Poisson bracket in classical field 
theory has for a long time been an open problem - as testified by the extensive literature on 
“multisymplectic Poisson brackets”, together with the fact that all these proposals suffer from 
serious defects. On the other hand, the functional approach does provide a good candidate 
which has come to be known as the Peierls - De Witt bracket and whose construction in a 
geometrical setting is now well understood. Here, we show how the basic “multisymplectic 
Poisson bracket” already proposed in the 1970s can be derived from the Peierls - De Witt 
bracket, applied to a special class of functionals. This relation allows to trace back most 
(if not all) of the problems encountered in the past to ambiguities (the relation between 
differential forms on multiphase space and the functionals they define is not one-to-one) and 
also to the fact that this class of functionals does not form a Poisson subalgebra. 
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1 Introduction 


The quest for a fully covariant Hamiltonian formulation of classical field theory has a long history 
which can be traced back to the work of Caratheodory [3], De Bonder [7] and Weyl [31] on the 
calculus of variations. From a modern point of view, one of the main motivations is the issue of 
quantization which, in traditional versions like canonical quantization as well as more recent ones 
such as deformation quantization, starts by bringing the classical theory into Hamiltonian form. 
In the context of mechanics, where one is dealing with systems with a finite number of degrees 
of freedom, this has led mathematicians to develop entire new areas of differential geometry, 
namely symplectic geometry and then Poisson geometry, whereas physicists have been motivated 
to embark on a more profound analysis of basic physical concepts such as those of states and 
observables. In the context of (relativistic) field theory, however, this is not sufficient since, 
besides facing the formidable mathematical problem of handling systems with an infinite number 
of degrees of freedom, we have to cope with new physical principles, most notably those of 
covariance and of locality. The principle of covariance states that meaningful laws of physics 
do not depend on the choice of (local) coordinates in space-time employed in their formulation: 
extending the axiom of Lorentz invariance in special relativity, it is one of the cornerstones of 
general relativity and underlies the modern geometrical approach to field theory as a whole. 
Equally important is the principle of locality, stating that events (including measurements of 
physical quantities) localized in regions of space-time that are spacelike separated cannot exert 
any influence on each other. Clearly, a mathematically and physically correct hamiltonian 
formalism for classical field theory should respect these principles: it should be manifestly 
covariant and local, as is the modern algebraic approach to quantum field theory; see, e.g.. 
Ref. |2]. 

As an example of a method that does not meet these requirements, we may quote the stan¬ 
dard hamiltonian formulation of classical field theory, based on a functional formalism in terms 
of Cauchy data: there, the mere necessity of choosing some Cauchy surface spoils covariance 
from the very beginning! To avoid that, a different approach is needed. 

Over the last few decades, attempts to construct such a different approach have produced a 
variety of proposals that, roughly speaking, can be assembled into two groups. 

One of these extends the geometrical tools that were so sucessful in mechanics to the situation 
encountered in field theory by treating spatial derivatives of fields on the same footing as time 
derivatives: in the context of a first order formalism, as in mechanics, this requires associating 
to each held component, say (/?*, not just one canonically conjugate momentum vr^ = dL/d^p^ 
but rather n canonically conjugate momenta = dL/d , where n is the dimension of 
space-time. (In mechanics, time is the only independent variable, so n = 1.) Identifying the 
appropriate geometrical context has led to the introduction of new geometrical entities now 
commonly referred to as “multisymplectic” and/or “polysymplectic” structures, and although 
their correct mathematical dehnition has only recently been completely elucidated m , the entire 
circle of ideas surrounding them is already reasonably well established, forming a new area of 
differential geometry; see mieHlIllSMI] for early references. 

A different line of thought is centered around the concept of “covariant phase space” laisiES], 
dehned as the space of solutions of the equations of motion: using this space to substitute the 
corresponding space of Cauchy data eliminates the need to refer to a specific choice of Cauchy 
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surface and has the additional benefit of providing an embedding into the larger space of all field 
configurations, allowing us to classify statements as valid “off shell” (i.e., on the entire space 
of field configurations) or “on shell” (i.e., only on the subspace of solutions of the equations of 
motion). 

Each of the two methods, the multisymplectic formalism as well as the covariant functional 
formalism, has its own merits and its own drawbacks, and experience has shown that best results 
are obtained by appropriately combining them. 

As an example to demonstrate how useful the combination of these two approaches can 
become, we shall in the present paper discuss the problem of giving an appropriate definition 
of the Poisson bracket, or better, the covariant Poisson bracket. From the point of view of 
quantization, this is a question of fundamental importance, given the fact that the Poisson 
bracket is expected to be the classical limit of the commutator in quantum field theory. Moreover, 
quantum field theory provides compelling motivation for discussing this limit in a covariant 
setting, taking into account that the (non-covariant) equal-time Poisson brackets of the standard 
hamiltonian formulation of classical field theory would correspond, in the sense of a classical 
limit, to the (non-covariant) equal-time commutators of quantum field theory, which are known 
not to exist in interacting quantum field theories, due to Haag’s theorem. 

Unfortunately, in the context of the multisymplectic formalism, the status of covariant 
Poisson brackets is highly unsatisfactory. This may come as a bit of a surprise, given the 
beautiful and conceptually simple situation prevailing in mechanics, where the existence of a 
Poisson bracket on the algebra C°°{P) of smooth functions on a manifold P is equivalent to 
the statement that P is a Poisson manifold and, as such, qualifies as a candidate for the phase 
space of a classical hamiltonian system: for any such system, the algebra of observables is just 
the Poisson algebra C°°{P) itself or, possibly, an appropriate subalgebra thereof, and the space 
of pure states is just the Poisson manifold P itself. In particular, this is true in the special case 
when P is a symplectic manifold, with symplectic form uj, say, and where the Poisson bracket of 
two functions f,g£ C°^{P) is the function {f,g} G C°°{P) defined by 

{f,g} = = w(X/,Ag) , (1) 

where Xj G X(P) denotes the hamiltonian vector field associated with / G C°°{P), i.e., 

ixf io = df . (2) 

This situation changes considerably, and for the worse, when we pass to the multisymplectic 
setting, where u is no longer a 2-form but rather an (n-|-l)-form and the hamiltonian vector 
field Xf is no longer associated with a function / but rather with an (n —l)-form /, n being 
the dimension of space-time. It can then be shown that equation ([2|) imposes restrictions not 
only on the type of vector field that is allowed on its Ihs but also on the type of differential 
form that is allowed on its rhs. Indeed, the validity of an equation of the form = df 

implies that the vector field X must be locally hamiltonian, i.e., we have Lxoj = 0, but it 

also implies that the form / must be hamiltonian, which by definition means that its exterior 
derivative df must vanish on all multivectors of degree n whose contraction with the (n -|- 1)- 
form uj is zero, and this is a non-trivial condition as soon as n > 1. (It is trivial for n = 1 

since u is assumed to be non-degenerate.) Thus it is only on the space of hamiltonian 

(n — l)-forms that equation ([T|) provides a reasonable candidate for a Poisson bracket. (Note, 
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however, that as a Lie algebra with respect to such a bracket, would have a huge 

center, containing the entire space of closed (n—l)-forms on P, since the linear map 

from (P) to Xlh{P) that takes f to Xf is far from being one-to-one: its kernel is precisely 
2 'ri-i(p) ) Anyway, the argument suggests that the transition from mechanics to field theory 
should somehow involve a replacement of functions by differential forms of degree n—1 - which 
is not completely unreasonable when we consider the fact that, in field theory, conservation laws 
are formulated in terms of conserved currents, which are closed (n —l)-forms. 

Unfortunately, this replacement leads to a whole bunch of serious problems, some of which 
are insurmountable. First and foremost, there is no reasonable candidate for an associative 
product on the space which would provide even a starting point for defining a Poisson 

algebra. Second, as has been observed repeatedly in the literature [T61I22H27] . the condition of 
being a locally hamiltonian vector field or a hamiltonian (n —l)-form forces these objects to 
depend at most linearly on the multimomentum variables, and moreover we can easily think of 
observables that are associated to forms of other degree (such as a scalar field, corresponding to 
a 0-form, or the electromagnetic field strength tensor, corresponding to a 2-form): this by itself 
provides enough evidence to conclude that hamiltonian (n —l)-forms constitute an extremely 
restricted class of observables and that setting up an adequate framework for general observables 
will require going beyond this domain. And finally, as has already also been noted long ago 
[ISllITlEolETHSg], the multisymplectic Poisson bracket defined by equation ([1]) fails to satisfy 
the Jacobi identity. In the case of an exact multisymplectic form (cj = —d9), this last problem 
can be cured by modifying the defining equation ([1]) through the addition of an exact (hence 
closed) term, as follows [12]: 

{f,9} = + d{ixj - ixf9 - ^Xg^Xj^) ■ (3) 

However, this does not settle any of the other two issues, namely 

• the lack of an associative product to construct a Poisson algebra; 

• the restriction to hamiltonian forms and forms of degree n—1, which leads to unreasonable 
constraints on the observables that are allowed, excluding some that appear naturally in 
physicists’ calculations. 

It should be mentioned here that these are long-standing problems: they have been recognized 
since the early stages of development of the subject (see, e.g., [T61I27] and also [251126] ') but have 
so far remained unsolved. 

A simple idea in this direction that has already been exploited is based on the observation 
that differential forms do admit a natural associative product, namely the wedge product, so one 
may ask what happens if, in the above construction, vector fields are replaced by multivector 
fields and (n—l)-forms by forms of arbitrary degree. As it turns out, this leads to a modified 
super-Poisson bracket, defined by a formula analogous to equation © Haul]. But it does not 
help to overcome either of the aforementioned other two issues. 

On the other hand, in the context of the covariant functional formalism, there is an obvious 
associative and commutative product, namely just the pointwise product of functionals, and 
apart from that, there also exists a natural and completely general definition of a covariant 
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Poisson bracket such that, when both are taken together, all the properties required of a Poisson 
algebra are satisfied; this bracket is known as the Peierls - De Witt bracket [SHinillSlES]. 

Thus the question arises as to what might be the relation, if any, between the covariant 
functional Poisson bracket, or Peierls - De Witt bracket, and the various candidates for multi- 
symplectic Poisson brackets that have been discussed in the literature, among them the ones 
written down in equations ([H) and (l3|). That is the question we shall address in this paper. 

In the remainder of this introduction, we want to briefly sketch the answer proposed here: 
details will be filled in later on. Starting out from the paradigm that, mathematically, classical 
fields are to be described by sections of fiber bundles, suppose we are given a fiber bundle P 
over a base manifold M, where M represents space-time, with projection p : P —M, and 
suppose that the classical fields appearing in the field theoretical model under study are sections 
(j) : M — > P of P (i.e., maps ip from M to P satisfying p o cp = idju), subject to appropriate 
regularity conditions: for the sake of definiteness, we shall assume here that all manifolds and 
bundles are “regular” in the sense of being smooth, while the regularity of sections may vary 
between smooth {C°°) and distributional To fix terminology, we define, for any section 

/ of any vector bundle V over P, its base support or space-time support, denoted here by supp/, 
to be the closure of the set of points in M such that the restriction of / to the corresponding 
fibers of P does not vanish identically, i.e.H 

supp/ = {x e M I f\p^ 7 ^ 0} . (4) 

Now suppose that / is a differential form on P of degree p, say, and that S is a closed p- 
dimensional submanifold of M, possibly with boundary, subject to the restriction that T, and 
supp/ should have compact intersectiorU so as to guarantee that the following integral is well 
defined, providing a functional j on the space of sections of P, 

(5) 

where p*f is of course the pull-back of / to M via p. Regarding boundary conditions, we shall 
usually require that if S has a boundary, it should not intersect the base support of /: 

dH n supp/ = 0 . (6) 

This simple construction provides an especially interesting class of functionals for various reasons, 
the most important of them being the fact that they are local, since they are simple integrals, over 
regions or more general submanifolds of space-time, of local densities such as, e.g., polynomials 
of the basic fields and their derivatives, up to a certain order 0 This is an intuitive notion 

no specification is given, it is tacitly assumed that we are dealing with smooth sections. 

^Using the abbreviation “supp” for the base support rather than the ordinary support (which would be a 
subset of P) constitutes a certain abuse of language, but will do no harm since the ordinary support will play no 
role in this paper. 

®Of course, this restriction is automatically satisfied if E is compact and also if / has compact base support. 
Moreover, given an arbitrary differential form / on P, we can always construct one with compact base support 
by multiplying with a “cutoff function”, i.e., the pull-back to P of a function of compact support on M. 

^To incorporate derivatives up to order r, say, of fields that are sections of some fiber bundle E over M, one 
has to define P using the r-th order jet bundle PE of E. 
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of locality for functionals of classical fields, but as has been shown recently, it can also be 
formulated in mathematically rigorous terms pQ. Either way, it is clear that the product of two 
local functionals of the form Q is no longer a local functional of the same form: rather, we get 
a “bilocal” functional associated with a submanifold of M x M and a differential form on P x P. 
Therefore, a mathematically interesting object to study would be the algebra of “multilocal” 
functionals which is generated by the local ones, much in the same way as, on an ordinary vector 
space, the algebra of polynomials is generated by the monomials. 

But the point of main interest for our work appears when we assume P to be a multi- 
symplectic fiber bundle m and M to be a Lorentz manifold, usually satisfying some additional 
hypotheses regarding its causal structure: more specihcally, we shall assume M to be globally 
hyperbolic since this is the property that allows us to speak of Cauchy surfaces. In fact, as is 
now well known, M will in this case admit a foliation by Cauchy surfaces dehned as the level 
sets of some smooth time function. However, it is often convenient not to fix any metric on M 
“a priori” since, in the context of general relativity, the space-time metric itself is a dynamical 
entity and not a fixed background held. Within this context, and for the case of a regular hrst- 
order hamiltonian system where helds are sections of a given conhguration bundle E over M 
and the dynamics is obtained from a regular hrst-order lagrangian via Legendre transform, it 
has been shown in [15] that one may use that structure to dehne the Peierls - De Witt bracket 
as a functional Poisson bracket on covariant phase space. Here, we want to show how, in the 
same context, multisymplectic Poisson brackets between forms, such as in equations m and dSj), 
can be derived from the Peierls - De Witt bracket between the corresponding functionals. For the 
sake of simplicity, this will be done for the case of (n—l)-forms, but we expect similar arguments 
to work in any degree. 

Concretely, we shall prove that given a hxed hypersurface S in M (typically, a Cauchy 
surface) and two hamiltonian (n—l)-forms / and g, we have 

> ^^,g} = ^T.,{f,g} > (7) 

where the bracket on the Ihs is the Peierls - De Witt bracket of functionals and the bracket {f,g} 
that appears on the rhs is a “multisymplectic pseudo-bracket” or “multisymplectic bracket” given 
by a formula analogous to equation © or to equation dSj). A more detailed explanation of this 
result will be deferred to the main body of the paper - last but not least because the construction 
requires the systematic use of both types of multiphase space that appear in field theory and that 
we refer to as “ordinary multiphase space” and “extended multiphase space”, respectively: they 
differ in that the latter is a one-dimensional extension of the former, obtained by including an 
additional scalar “energy type” variable. Geometrically, extended multiphase space is an affine 
line bundle over ordinary multiphase space, and the hamiltonian P of any theory with this type 
of “field content” is a section of this affine line bundle. Moreover, each of these two multiphase 
spaces comes equipped with a multisymplectic structure which is exact (i.e., the multisymplectic 
form is, up to a sign introduced merely for convenience, the exterior derivative of a multicanonical 
form), naturally dehned as follows. First, one constructs the multisymplectic form ui and the 
multicanonical form 6 on the extended multiphase space by means of a procedure that can be 
thought of as a generalization of the construction of the symplectic structure on the cotangent 
bundle of an arbitrary manifold. Then, the corresponding forms on the ordinary multiphase 
space are obtained from the previous ones by pull-back via the hamiltonian P: therefore, they 
will in what follows be denoted by and by to indicate their dependence on the choice 
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of hamiltonian. We can express this by saying that the multisymplectic structure on extended 
multiphase space is “kinematical” whereas that on ordinary multiphase space is “dynamical”. 
Correspondingly, we shall refer to the brackets defined by equations & and m on extended 
multiphase space as “kinematical multisymplectic brackets” and to the brackets defined by the 
analogous equations 

{/) 9 } = ixg'^Xf^H — (^) 

with 

ixf<^H = df (9) 

and 

{/)9} = 'i'Xg'^Xf^n + d(^ixgf — ixf9 — ixgixfdnj (^O) 

as “dynamical multisymplectic brackets”. In both cases, the brackets defined by the simpler 
formulas ([T|) and ([SD are really only “pseudo-brackets” because they fail to satisfy the Jacobi 
identity, and the correction terms that appear in equations ([3]) and (fTOt) are introduced to cure 
this defect. At any rate, what appears on the rhs of equation ([7|) above is the dynamical bracket 
on ordinary multiphase space and not the kinematical bracket on extended multiphase space - 
in accordance with the fact that the Peierls - De Witt bracket itself is dynamical. 

We conclude this introduction with a few comments about the organization of the paper. 
In Section 2, we set up the geometric context for the functional calculus in classical field theory, 
introduce the class of local functionals to be investigated and give an explicit formula for their 
first functional derivative. In Section 3, we present a few elementary concepts from multisym¬ 
plectic geometry, which is the adequate mathematical background for the covariant hamiltonian 
formulation of classical field theory. In Section 4, we combine the two previous sections to 
formulate, in this context, the variational principle that provides the dynamics and derive not 
only the equations of motion (De Bonder--Weyl equations) but also their linearization around 
a given solution (linearized De Bonder - Weyl equations), with emphasis on a correct treatment 
of boundary conditions. In Section 5, we present the classification of locally hamiltonian and 
exact hamiltonian vector fields on multiphase space. Section 6 contains the main result of the 
paper on the connection between multisymplectic Poisson brackets and the functional Poisson 
bracket of Peierls and DeWitt. Finally, Section 7 provides further discussion of this result, its 
implications and perspectives for future investigations. 

The paper presents a substantially revised and expanded version of the main results contained 
in the PhD thesis of the second author [33], which was elaborated under the supervision of the 
hrst author. 


2 Geometric setup for the functional calculus 

We begin by collecting some concepts and notations that we use throughout the article. 
As already mentioned in the introduction, classical fields are sections of fiber bundles over space- 
time, so our starting point will be to fix a fiber bundle P over the space-time manifold M (not 
necessarily endowed with a fixed metric, as mentioned before), with projection p : P — M. 
The space of held conhgurations C is then the space of (smooth) sections of P, or an appropriate 
subspace thereof, 

e C r°°(P) , (11) 
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whose elements will, typically, be denoted by cj). Formally, we can view this space as a mani¬ 
fold which, at each point 4>, has a tangent space T^Q and, similarly, a cotangent space T^Q- 
Explicitly, denoting by the vertical bundle of P, pulled back to M via (j), 

= riYeiP) , ( 12 ) 

and by V® its twisted dual, defined by taking the tensor product of its ordinary dual with the 
line bundle of volume forms over the base space, 

V® = v; ® A^T*M , (13) 

we have that, according to the principles of the variational calculus, T^C is the space of smooth 
sections of V^, or an appropriate subspace thereof, 

c r-(E^), (14) 

whose elements will, typically, be denoted by 64> and called variations of (j), whereas T^C is the 
space of distributional sections of V®, or an appropriate subspace thereof, 

r;e c r—(F®). (i5) 

The reader will note that in equations (jlip . (I14p and ()15p . we have required only inclusion, rather 
than equality. One reason is that the system may be subject to constraints on the fields which 
cannot be reduced to the simple statement that they should take values in some appropriate 
subbundle of the original bundle (this case could be handled by simply changing the choice of 
the bundle P). But even for unconstrained systems, which are the only ones that we shall be 
dealing with in this paper, there is another reason, namely that we have not yet specified the 
support properties to be employed. One obvious possibility is to set 

r^e = r-(F^) , T;e = r—( f®) , (le) 

which amounts to allowing only variations with compact support. At the other extreme, we may 
set 

T^e = r-(F^) , T;e = r--(F®). (i7) 

And finally, there is a third option, specifically adapted to the situation where the base space is 

a globally hyperbolic lorentzian manifold and adopted in m, which is to take 

= r-(F^) , r;e = r,--(F®), (is) 

where the symbols “sc” and “tc” indicate that the sections are required to have spatially compact 
support and temporally compact support, respectively. These options correspond to different 
choices for the support properties of the functionals that will be allowed. 

Generally speaking, given a functional F on C, we define its base support or space-time 
support, denoted here simply by supp F, as follows [T]: 

There exists an open neighborhood Ux of x in M 
X ^ supp F such that for any two field conhgurations (j)i,4>2 F Q (19) 

satisfying (j)^ = (p 2 on M \Ux , F[(j)^] = F[4>2\- 
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This definition implies that supp F is a closed subset of M since its complement is open: it is 
the largest open subset of M such that, intuitively speaking, F is insensitive to variations of 
its argument localized within that open subset. It also implies that the functional derivative of F 
(if it exists) satisfies 

= 0 if supp F n supp 5i?!) = 0. (20) 

For later use, we note that the functional derivative will often be expressed in terms of a (formal) 
variational derivative: 

= j ^[4>\{x) 54>\x) . (21) 

Jm ^(P 

Typically, as always in distribution theory, the functional derivative will be well defined on 
variations dcp such that supp F n suppdt/) is compact. Thus if no restrictions on the space-time 
support of F are imposed, we must adopt the choice made in equation (llOp . At the other extreme, 
if the space-time support of F is supposed to be compact, we may adopt the choice made in 
equation m- And finally, the choice made in equation (fT8]l is the adequate one for dealing 
with functionals that have temporally compact support, i.e., space-time support contained in the 
inverse image of a bounded interval in M under some global time function: the typical example 
is that of a local functional of the form given by equation ([5]) when S is some Cauchy surface. 
More generally, note that for local functionals of the form given by equation ([5]), we have 

suppJ^j^j = S n supp/ . (22) 

However, it should not be left unnoticed that the equality in equations (fT6]) - (fT8]l and, possibly, 
in equation (1221) . can only be guaranteed to hold for non-degenerate systems, since in the case 
of degenerate systems, there will be additional constraints implying that we must return to the 
option of replacing equalities by inclusions, as before. 

In what follows, we shall make extensive use of the fact that variations of sections can always 
be witten as compositions with projectable vector fields, or even with vertical vector fields, on 
the total space P. To explain this, recall that a vector field X on the total space of a fiber 
bundle is called projectable if the tangent map to the bundle projection takes the values of X 
at any two points in the same fiber to the same tangent vector at the corresponding base point, 
i.e., 

■ ^{Pi) = Tp^P ■ ^{Pi) for e P such that p{p^) = p{p 2 ) ■ (23) 

This is equivalent to requiring that there exists a vector field Xj^ on the base space which is 
yO-related to X, 

^m(^) = '^pP ' ^(p) for p G P such that p{p) = m . (24) 

In particular, X is called vertical if Xm vanishes. Now note that given any projectable vector 
field X on P, we obtain a functional vector field X on C whose value at each point /> G C is the 
functional tangent vector X[/>] G T^Q, denoted in what follows by defined as 

SxcP = X{<P) - TcfiX^) , (25) 


or more explicitly. 


6x4>{m) = X{(j){m)) — T^(f{Xxi{m)) for m^M . 


(26) 


Conversely, it can be shown that every functional tangent vector can be obtained in this way 
from a vertical vector field X on P, i.e., given a section 64> of (j)*(VerP), there exists a vertical 
vector field X on P representing it in the sense that 6(p is equal to Sx4>- (To do so, we can 
apply the implicit function theorem to construct, for any point m in M, a system of local 
coordinates (x^,?/") for P around (j)[m) in which p corresponds to the projection onto the 
first factor, (x^, ?/“) i—>■ and (p corresponds to the embedding in the first factor, x^ i—>■ (x^, 0). 
Moreover, in these coordinates, 64> is given by functions 50“ (x^) whereas X is given by functions 
X“(x/^,/), so we may simply define an extension of the former to the latter by requiring the 
X“ to be independent of the , setting X“(x^,?/^) = 50“(x^) in a neighborhood of the origin 
in y-space and then using a smooth cutoff function in y-space.) 

Of course, the reader may wonder why, in this context, we bother to allow for projectable 
vector fields rather than just vertical ones. The point is that although vertical vector fields 
are entirely sufficient to represent variations of sections, we shall often encounter the con¬ 
verse situation in which we consider variations of sections induced by vector fields which are 
not vertical but only projectable, such as the hamiltonian vector fields appearing in equations 

(PP)-® and (I5])-([ini). 

Regarding notation, we shall often think of a projectable vector field as a pair X = {Xp, Xm) 
consisting of a vector field Xp on the total space P and a vector field Xm on the base space M, 
related to one another by the bundle projection: then equations (|25|) and (|26p should be written 
as 

= Xp{P) -TP{Xm) , (27) 

and 

6x<P{m) = Xp{4>{m)) — T^cj){XM{m)) for m^M , (28) 

respectively. The same argument as in the previous paragraph can then be used to prove the 
following 

Lemma 1 Let p he a section of a fiber bundle P over a base manifold M. Given any vector 
field Xj^ on M, there exists a projectable vector field xf, on P which is fi-related to X^, i.e., 
satisfies Xp(0) = Tcj){XM), and then we have (pi*{i^^a) = ix^iP*^), for any differential form 
a on P. 

As an example of how useful the representation of variations of sections of a fiber bundle 
by composition with vertical vector fields or even projectable vector fields can be, we present 
explicit formulas for the first and second functional derivative of a local functional of the type 
considered above - beginning with a more detailed definition of this class of functionals. 

Definition 1 Given a fiber bundle P over an n-dimensional base manifold M, let he a 
p-dimensional submanifold of M, possibly with boundary dH, and f he a p-form on the total 
space P such that the intersection of S with the base support of f is compact. The local 
functional associated to S and f is the functional j : C —>■ M on the space C C T°°{P) 
of field configurations defined by 

Xj:,f[P] = forPGQ. (29) 
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These functionals are differentiable, and their derivative is given by a completely explicit formula: 


Proposition 1 Given a fiber bundle P over an n-dimensional base manifold M, let T, be a 
p-dimensional submanifold of M, possibly with boundary dTi, and f be a p-form on the total 
space P such that the intersection of S with the base support of f is compact. Then the local 
functional j associated to S and f is differentiable, and representing variations of sections 
of P in the form 6x4> where X = {Xp,X^) is a projectable vector field, its functional derivative 
is given by the formula 

^ {^*{Lxj) - Lx^{^*f)) for , (30) 

where Lz denotes the Lie derivative along the vector field Z. 

Remark 1 Under the boundary condition that the intersection of with the base support 
of / is empty, equation (IHU]) can be rewritten as follows: 

Tfj\f\-bxf> = J^{f*{ix,df) -ix^{f*df)) for (31) 

The same equation holds when this boundary condition is replaced by the requirement that 6^4> 
should vanish on dT,. 


Proof: Recall first that for any functional F on 6, its functional derivative at <() € C along 

bcj) €T^Q is defined by 




A=0 


where the a Q constitute a smooth one-parameter family of sections of P such that 


= 


A=o dX 


A=0 


Fixing (p and dp and choosing a projectable vector field X = {Xp,Xj^) that represents 5p as 
bxP-, according to equation (1251) . consider its flow, which is a (local) one-parameter group of 
(local) automorphisms = (<l>p^,d>^^) such that 


Xp = 


d ^ 


A=0 


and X]^ = 


d\ 


$ 


M,A 


A=0 


This allows us to take the one-parameter family of sections of P to be given by the one- 
parameter group of automorphisms <1>^, according to 


^M,X ■ 
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A=0 


Now we are ready to calculate: 



To derive equation (lai) from equation dSQ]), it suffices to apply standard formulas such as 
Lz = diz + iz d and the fact that d commutes with pull-backs, together with Stokes’ theorem, 
and use the boundary condition ([6]) to kill the resulting two integrals over dT,. The same 
argument works when is supposed to vanish on OS, since we may then arrange Xj^ to 
vanish on OS and Xp to vanish on □ 


3 Multiphase spaces and multisymplectic structure 

As has already been stated before, the bundle P appearing in the previous two sections, rep¬ 
resenting the multiphase space of the system under consideration, will be required to carry 
additional structure, namely a multisymplectic form. There has been much debate and even 
some confusion in the literature on what should be the “right” definition of the concept of a 
multisymplectic structure, but all proposals made so far can be subsumed under the following 

Definition 2 A multisymplectic form on a manifold P is a differential form ui on P of 
degree n-|-l, say, which is closed, 

du! = 0 , (32) 

and satisfies certain additional algebraic constraints - among them that of being non-degenerate, 
in the sense that for any vector field X on P, we have 

ixuj = 0 X = 0 . (33) 

Of course, this definition is somewhat vague since it leaves open what other algebraic constraints 
should be imposed besides non-degeneracy. One rather natural criterion is that they should be 
sufficient to guarantee the validity of a Darboux type theorem. Clearly, when n = l, the above 
definition reduces to that of a symplectic form, and no additional constraints are needed. But 
when n > 1, which is the case pertaining to field theory rather than to mechanics, this is no 
longer so. An important aspect here is that P is not simply a manifold but rather the total 
space of a fiber bundle over the space-time manifold M, which is supposed to be n-dimensional, 
so one restriction is that the degree of the form u is linked to the space-time dimension. Another 
restriction is that cu should be (n —l)-horizontal, which means that its contraction with three 
vertical vector fields vanishes: 


ixiy'i-z^ — 0 foi" vertical . 


(34) 
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And finally, there is a restriction that, roughly speaking, guarantees existence of a “sufficiently 
high-dimensional” lagrangian subbundle of the tangent bundle of P, but since we shall not need 
it here, we omit the details: they can be found in Ref. m- 

The main advantage of the definition of a multisymplectic form as given above is that we 
can proceed to discuss a number of concepts which do not depend on the precise nature of the 
additional algebraic constraints to be imposed. For example, a vector field A" on P is said to be 
locally hamiltonian if is closed, that is, if0 


LxUJ = 0 , (35) 

and is said to be globally hamiltonian or simply hamiltonian if i^u} is exact, that is, if there 
exists an (n—l)-form / on P such that 


ixU} = df . (36) 

Reciprocally, an (n—l)-form / on P is said to be hamiltonian if there exists a vector held X 
on P such that equation (I36h holds: this condition is trivially satished when n = 1 but not 
when n > 1. Note that due to non-degeneracy of w, X is uniquely determined by /, and will 
therefore often denoted hy Xf, whereas / is determined by X only up to addition of a closed 
form: despite this (partial) ambiguity, we shall say that X is associated with / and / is associated 
with X. In the special case when w is exact, i.e., we have 

io = -de , (37) 

where 9 is an appropriate n-form on P called the multicanonical form, a vector held X on P 
is said to be exact hamiltonian if 

Lx9 = 0 . (38) 

In this case, of course, the associated hamiltonian form can be simply chosen to be ix^, since 
dix9 = LxO — ixdO = ix^- In particular, this happens when P is the total space of a vector 
bundle (over some base space E, say, which in turn will be the total space of a hber bundle over 
the space-time manifold M), provided that w is homogeneous of degree one with respect to the 
corresponding Euler vector field or scaling vector field S, i.e., 

L-^uj = uj , (39) 

since we may then dehne 9 by 

9 = — ij^uj . (40) 

Moreover, we can then employ S to decompose vector helds and differential forms on P ac¬ 
cording to their scaling degree, and as we shall see below, this turns out to be extremely useful 
for the classihcation of hamiltonian vector helds (whether locally or globally or exact) and of 
hamiltonian forms. 

The standard example of this kind of structure is provided by any (hrst order) hamiltonian 
system obtained from a (hrst order) lagrangian system via a non-degenerate covariant Legendre 
transformation. In this approach, one starts out from another hber bundle over M, denoted 
here by E and called the configuration bundle: the relation between the two bundles E and P 

^Throughout this paper, we shall make use of Cartan’s formula Lx — dix + ixd without further mention. 
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is then established by taking appropriate duals of hrst order jet bundles. Namely, consider 
the hrst order jet bundle of E, denoted simply by JE, which is both a hber bundle over M 
(with respect to the source projection) and an affine bundle over E (with respect to the target 
projection), together with its difference vector bundle, called the linearized hrst order jet bundle 
of E and denoted here by JE, which is both a hber bundle over M (with respect to the source 
projection) and a vector bundle over E (with respect to the target projection), and introduce 
the corresponding duals: the affine dual J®E of JE and the usual linear dual J®E of JE. 
In what follows, we shall refer to the latter as the ordinary multiphase space and to the former 
as the extended multiphase space of the theory. As it turns out and has been emphasized since 
the beginning of the “modern phase” of the development of the subject in the early 1990s 
(see, e.g., 0), both of these play an important role since not only are both of them fiber 
bundles over M and vector bundles over E, but J®E is also an affine line bundle over J®E, 
and the dynamics of the theory is given by the choice of a hamiltonian, which is a section 
H : J®E —>■ J®E of this affine line bundle. Moreover, and this is of central importance, 
both of these multiphase spaces carry a multisymplectic structure. Namely, J®E comes with a 
naturally defined multisymplectic form of degree n+1, denoted here by uj, which (up to sign) is 
the exterior derivative of an equally naturally dehned multicanonical form of degree n, denoted 
here hy 9 : ui = —dO, and if we choose a hamiltonian Ji : J®E —^ J®E, we can pull them back 
to obtain a corresponding multicanonical form 

9^ = 'H*9 (41) 

and a corresponding multisymplectic form 

= H*uj (42) 

on J®E: again, ujn = —ddy^. Thus the main difference between the extended and the ordinary 
multiphase space is that 9 and u are universal and “cinematical”, whereas 9-^ and luxl are 
“dynamical”. In terms of local Darboux coordinates {x^,q^,p^) for J®E and 
for J®E induced by the choice of local coordinates for M, g* for the typical fiber Q oi E and 
a local trivialization of E over M, we have 



9 

= pf dq^ A d'^x^ + p d^x , 

(43) 

and 

OJ = 

dg* A dp^ A d^'x^ — dp A d^x , 

(44) 

so that writing Ti = 

-Hd^x, 





= pf dq^ A d^'Xfj, — H d^x , 

(45) 

and 


dq^ A dp^ A d'^Xfj, + dH a d'^x , 

(46) 

or more explicitly 





< 

< 

II 

d'^x^ + ^ dg* A d^x + dpf a d”x . 

ag* opf 

(47) 
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where is the (local) (n—l)-form obtained by contracting the (local) volume form d^x with 
the local vector field = d/dx^] for more details, including a global definition of 0 that does 
not depend on any of these choices, we refer to 

4 Variational principle and equations of motion 

The fundamental link that merges the functional and multisymplectic formalisms discussed in 
the previous two sections into one common picture becomes apparent when the construction of 
functionals of fields from forms on multiphase space outlined in the introduction is applied to the 
multicanonical re-form on (ordinary) multiphase space: this provides the action functional S 
of the theory, defining the variational principle whose stationary points are the solutions of the 
equations of motion. Indeed, the action functional is really an entire family of functionals 
on the space 6 of field configurations f (see equation (fTT]) i. given by 

= [ r Oh, (48) 

JK 

where K runs through the compact submanifolds of M which are the closure of their interior 
in M and have smooth boundary dK. 

Within this setup, a section <() in S is said to be a stationary point of the action if, for any 
compact submanifold K oi M which is the closure of its interior in M and has smooth boundary 
dK, (p becomes a critical point of the functional Sr restricted to the (formal) submanifold 

G C I ^\dK = 4>\dK} (49) 

of C, or equivalently, if the functional derivative Snip] of at p vanishes on the subspace 

eT^e\6p = 0 on dK} (50) 

of T^Q- As is well known, this is the case if and only if p satisfies the corresponding equations 
of motion, which in the present case are the De Bonder - Weyl equations; see, e.g., [T^ITSIITU] . 
Globally, these can be cast in the form 


p*{ixOJ^ — 0 

for any vertical vector field A on P ’ 


(51) 


or even 


= 0 

for any projectable vector field X = {Xp,Xj^) 
whereas, when written in terms of local Darboux coordinates (x^, q^,pf) as before, they read 


(52) 


n * dH 


(53) 


®It should be noted that whereas the form uj is always non-degenerate, the form ujpi is degenerate in mechanics 
[n = 1) and non-degenerate in field theory (n > 1). 
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where P = J®E, cp = {(p,7r) and P = —Hd^x. Similarly, given such a solution (/>, we shall 
say that a section 64> in T^Q is an infinitesimal stationary point of the action if it is formally 
tangent to the (formal) submanifold of solutions, or equivalently, if it satisfies the corresponding 
linearized equations of motion, which in the present case are the linearized De Bonder ~ Weyl 
equations. Globally, representing Scj) in the form 5xfi where X = {Xp,Xj^) is a projectable 
vector field, these can be cast in the form 


(j)* {iyLxp^Xy) — 0 

for any vertical vector field B on P 


(54) 


or even 


for any projectable vector field Y = (Yp,Yj^) 
whereas, when written in terms of local Darboux coordinates (x^, q^,Pi) as before, they read 


(55) 


^ ^ d^H , , ^ ^ 


dp^j dpr 




(56) 


dq^ dq'- 


dpj dq^ 


where P = J®E, cf = (v^, vr), 5(j) = {5ip, 5 tt) and P = —H d^x. 

Proof: For the first part (concerning the derivation of the full equations of motion from the 

variational principle), we begin by specializing equation (I.SOD to vertical vector fields X on P 
(i.e., setting Xj^ = 0 and replacing Xp hy X), with f = and using standard facts such as 
the formula Lz = diz + iz d or that d commutes with Lie derivatives and pull-backs, together 
with Stokes’ theorem, to obtain that, for any vertical vector held X on P, 


' dx4> — — [ fi*{ix^n) + [ 4’*{'^xd'hL) ■ 

■JK JdK 

Obviously, condition (j5ip implies that this expression will be equal to zero for all vertical vector 
helds X on P which vanish on P\dK- Conversely, it follows from Lemma [2] below that if this 
is the case, then condition (|51ll holds. Moreover, it is easily seen that this condition is really 
equivalent to the condition 


for any projectable vector held X = (Xp,Xj^^) 


but it so happens that the form is identically zero, for dimensional reasons. Finally, 

a simple calculation shows that equation ()5ip . when written out explicitly in local Darboux 
coordinates, can be reduced to the system ([53p . 

For the second part (concerning the linearization of the full equations of motion around a given 
solution ())), we proceed as in the proof of Proposition [TJ hxing and dcj), suppose we are given 
a smooth one-parameter family of sections (j)^ a C of P such that 


= </’A 


A=0 


S4> = ^<l>. 


A=0 
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as well as a projectable vector field X = {Xp,X^) that represents Scj) as according to 

equation (1251) . together with its flow, which is a (local) one-parameter group of (local) automor¬ 
phisms <1>;^ = ($p;^,<l>^^) such that 




A=0 


and Xj^ = 


d\ 




M,A 


A=0 


allowing us to take the one-parameter family of sections (j))^ of P to be given by the one-parameter 
group of automorphisms according to 


4 >\ — ‘^p,a' 




-1 

M,A • 


Then for any vertical vector held Y on P, we have 


d_ 

dX 


^h) 


A=0 


d 


> O $ 


-1 

M,A 


{iy 


A=0 

d 


A=0 


(t>*{LxpiY^H) - Lxu{4>*{iY^'H)) 

4>*{iY^Xp^'H) + ‘P* {'i'[Xp,Y]^'^) ~ d‘XM{'P*^^Y ^h)) ; 


where the last two terms vanish according to equation ()51l) . and the same argument holds with 
Y replaced by Yp where Y = {Yp,Yj^) is any projectable vector field. Finally, an elementary 
but somewhat lengthy calculation shows that equation (I54p . when written out explicitly in local 
Darboux coordinates, can be reduced to the system (1561) . provided 4> satisfies the system (f53]l . 

□ 


The lemma we have used in the course of the argument is the following. 


Lemma 2 Given a fiber bundle P over an n-dimensional base manifold M, let be a section 
of P and OL be an {n + l)-form on P such that, for any compact submanifold K of M which is 
the closure of its interior in M and has smooth boundary dK, and for any vertical vector field 
X on P that vanishes on P\qk together with all its derivatives, the integral 



vanishes. Then the form (p*{ixa) vanishes identically, for any vertical vector field X on P (not 
subject to any boundary conditions). 


Proof: Suppose X is any vertical vector field on P and m is a point in M where cjfpixOt) does 

not vanish. Choosing an appropriately oriented system of local coordinates in M around m, 
we may write 4>* (ijs^o) = a d^x where a is a function that is strictly positive at the coordinate 
origin (corresponding to the point m) and hence, for an appropriate choice of sufficiently small 
positive numbers 5 and e, will be ^ e on Bs (here we denote by and by B^ the open and 
closed ball of radius r around the coordinate origin, respectively). Choosing a test function x 
on M such that 0 ^ x ^ 1; suppx C Bs and x = ^ on Bs/ 2 , lifted to P by pull-back and 


16 








multiplied by X to give a new vertical vector field yX on P that vanishes on PIqb^ together 
with all its derivatives, we get 



d"'x xix)a{x) 





d'^x a(x) 


^ e vo 1 (Bs/ 2 ) > 0 , 


which is a contradiction. 


□ 


To summarize, the space of stationary points of the action, denoted here by S and considered 
as a (formal) submanifold of the space of held conhgurations C, can be described in several 
equivalent ways: we have 

§ = e I 0 is stationary point of the action} , (57) 

or 

§ = {(() € C I 0 satishes the equations of motion (f53]l } , (58) 

or 

S = {(/> G C I ~ projectable vector held X = (Xp,X^)} , (59) 

or 

S = {(() G C I (p*{ix ^n) — 0 vertical vector held X on P] . (60) 


This space § plays a central role in the functional approach: it is widely known under the name 
of covariant phase space. Moreover, given in S and allowing X = {Xp,Xj^) to run through 
the projectable vector helds, the (formal) tangent space T.S to S at (() is 


or 


or 


T^S = {dcp € T^Q I 6(j) is inhnitesimal stationary point of the action} , (61) 

r^S = G r^C I d(j) satishes the linearized equations of motion (f56]l } , (62) 


r^s = {6x<peT^e 


i'^Yp^Xp^H 


) = 0 


for any projectable vector held Y = (Yp, i 




(63) 


or 


r^s = {dxcPeT^e 


(j)*{iYL 


Xp^n. 


= 0 


for any vertical vector held Y on P 


}• 


(64) 


5 Locally hamiltonian vector fields 

The results of the previous section, in particular equations (IUH|1 and HMD , provide strong moti¬ 
vation for studying projectable vector helds X = {Xp,Xj^) on (ordinary) multiphase space 
which are locally hamiltonian (that is, such that Xp is locally hamiltonian), since they imply 
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that each of these provides a functional vector field X on covariant phase space defined by a 
simple algebraic composition rule: 

X[(l)] = 6x(t> for (/>€§. (65) 

As we shall see in the next section, this functional vector field is also hamiltonian with respect 
to the natural symplectic form on covariant phase space to be presented there. But before doing 
so, we want to address the problem of classifying the locally hamiltonian vector fields and, along 
the way, also the globally hamiltonian and exact hamiltonian vector fields on multiphase space. 
For the case of extended multiphase space, using the forms oj and 6, this classification has been 
available in the literature for some time naE], even for the general case of multivector fields. 
But what is relevant here is the corresponding result for the case of ordinary multiphase space, 
using the forms uj-^ and 9^, for a given hamiltonian T-L. As we shall see, there is one basic 
phenomenon common to both cases, namely that when re> 1 , the condition of a vector field to 
be locally hamiltonian imposes severe restrictions on the momentum dependence of its compo¬ 
nents, forcing them to be at most affine (linear -|- constant): this appears to be a characteristic 
feature distinguishing mechanics (n = l) from field theory (n>l) and has been noticed early by 
various authors (see, e.g., [I61I27]) and repeatedly rediscovered later on (see, e.g.. 

However, despite all similarities, some of the details depend on which of the two types of multi¬ 
phase space we are working with, so it seems worthwhile to give a full statement for both cases, 
for the sake of comparison. 

We begin with an explicit computation in local Darboux coordinates for J®E 

and (x^, q^,Pi,p) for J®E induced by the choice of local coordinates for M, 5 * for the typical 
fiber Q oi E and a local trivialization of E over M, supposing that we are given a hamiltonian 
Ti : J®E —> J®E which we represent in the form T-L = — Hd^x, as usual. Starting out from 
equations dH-dUl), we distinguish two cases: 


Extended multiphase space 

Using equations (|43p - (l44p and writing an arbitrary vector held X on J®E as 


X = 


f) -3 3 

^ I _ _L 

dxf^ dq^ * dp^ 




( 66 ) 


we hrst note that X will be projectable to E if and only if the coefficients X^ and X® do not 
depend on the energy variable p nor on the multimomentum variables and will be projectable 
to M if and only if the coefficients X^ do not depend on the energy variable p nor on the 
multimomentum variables nor on the position variables q^. Next, we compute 


= X'^ dq^ A dp^ A d^x^u — X^ dq'‘ a d'^x^ + X® dp^ a d'^x^ 

-\- X^ dp A d'^x^ — Xg d'^x , 


and 

Then we have 


ix9 = (pfX® + pX'®) d'^x^ - pfX'^ dq^Ad^x 




(67) 

( 68 ) 
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Theorem 1 Any locally hamiltonian vector field on J®E is projectable to E. Moreover, given 
any vector field X on J®E which is projectable to E, X is locally hamiltonian if and only if 


1. If N>1, X is also projectable to M, i.e., the coefficients X^ do not depend on the position 
variables q^. 

2. The coefficients X^ and Xq can be expressed in terms of the previous ones and of new 
coefficients Xf_ which, once again, do not depend on the energy variable p nor on the 
multimomentum variables pf, according to 


= -p 


dX>‘ „ dXi 


dq^ 


Pi 


dq^ 


+ Pi 


dX^^ 

dx’^ 


Pi 


dx^ dxt 

dx'^ dq^ 


X, = -p 


dXf^ „ dX^ dXt 


Pi 


+ 


dx>^ * ^x^^ ^x^^ 

Finally, X is exact hamiltonian if and only if, in addition, the coefficients Xf_ vanish. 


(69) 

(70) 


Proof: The proof is carried out by “brute force” computation and a term by term analysis 

of the coefficents that appear in the Lie derivative of w along X: it will be omitted since an 
explicit proof of a more general statement along these lines can be found in HU. □ 


Ordinary multiphase space 


Using equations dlUl-dlTl) and writing an arbitrary vector held X on J®E as 

3 -3 3 

^ — + X^ — + xi^ ^ , 

dx>^ dq^ dp^ 


(71) 


we hrst note that X will be projectable to E if and only if the coefficients X^ and X® do 
not depend on the multimomentum variables pf and will be projectable to M if and only if 
the coefficients X^ do not depend on the multimomentum variables p'f. nor on the position 
variables q^. Next, we compute 


ix^H = X'’' dq^ A dpf A d"'Xfj,y — + X^ 'gXj ^ 

+ + (^‘ I? + w) • 

and 

ixdn = {PiX^ - HX^) d^x^ - pfX" dq^Ad^Xf,, . 

Then we have 


(72) 


(73) 


Theorem 2 Any locally hamiltonian vector field on J®E is projectable to E, except possibly 
when N = 1 and n = 2, and any exact hamiltonian vector field on J®E is projectable to E. 
Moreover, given any vector field X on J®E which is projectable to E, X is locally hamiltonian 
if and only if 
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1. If N>\, X is also projectable to M, i.e., the coefficients do not depend on the position 
variables q^. 

2. The coefficients X^ can be expressed in terms of the previous ones and of new coefficients 
Xf_ which, once again, do not depend on the multimomentum variables p^, according to^ 


Xf = H 


^X^^ 


dX^ 


dq^ dq^ ' dx^ dx^ ' dq^ 

3. The coefficients X^, X^ and Xf satisfy the compatibility condition 


T- + Pi 


dXi^ 


dX^ 


Pi 


+ 


dX^^ 


(74) 


dH 


dH 


dH 


dxf^ ^ dq^^ ^ 


Xf = - H 


dXf^ 


+ Pi 


dx^ dxf 


dxf^ * dx>^ dx>^ 
Finally, X is exact hamiltonian if and only if, in addition, the coefficients Xf_ vanish. 


(75) 


Maybe somewhat surprisingly, the case where A^ = 1 and n = 2, concerning the theory of a single 
real scalar field in two space-time dimensions, is somewhat exceptional in that it allows for 
locally hamiltonian vector fields which fail to be projectable and are not covered by the above 
classification theorem; we shall address this question in Remark [2] below. 

Proof: As in the case of Theorem 1, the proof is carried out by “brute force” computation. 

First, we apply the exterior derivative to equation ([72]) and collect terms to get 


= (^(^r + + tK'i I 


dx 


+ <5; 


d 


^ ^dr,x^ - 

OPi 


V \ M' 




dX'^ dX° 


dx'^ dx^ 


dH 


dq^ 


d 


dq^ 


a \rk 


dH 




dpi 


d 


dq^ 

d 


dH 


- — Af + A dg* a d 


(9g® 


X, 


dpK\ A 


+ — ( d? A' - A'” ““ ' 


djj' 

— ) dpt A dpf A d'^x^ 


dpi 


dX^ 7 fc 1 l j A in 
+ ITT A dg A dpi a d x^^ 


dq^ 

dX^ 

dpi 


dq^ A dpi A dpi; A d'^xxo 


J . . an„ 


X„ 


^Note that the first term in equation dZH) is absent as soon as A^> 1. 
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Numbering the terms in this equation from 1 to 7, we begin by analyzing terms no. 7, 5 and 6 . 
Obviously, when X is projectable to E and the coefficients satisfy the condition stated in 
item 1 . of the theorem, these terms vanish identically (term no. 6 is absent when = 1 ), so 
what we need to analyze is the converse statement. 


Term No. 7: For any choice of indices and mutually different indices 

pi,, Pn- 2 , contracting L^ujy^ with the multivector field a a a a ... a 
gives the relation 

dX^ : dX^ 


5^ 


dPi 


At ^uapi...pn-2 


dp 


u ^^ap\...pn-2 


(76) 


Now if to begin with, we fix only the indices i and p, together with some other index p, 
we can always choose the remaining free indices in this equation to be such that m=j and 
i',pi,..., Pn -2 are all mutually different and ^ p: this reduces the Ihs to the expression 
EdXP!dp^, while the rhs vanishes if we take m^i, which is possible as soon as > 1 . 
Thus we conclude that 

dXP 


dp^ 


= 0 


except perhaps when N = l. But even when A^ = l, where equation ([76]l reduces to|^ 

dX^ dX^ 




p ^vapi...p„-2 


dp 


V PCrpx...pn-2 ’ 


(77) 


this conclusion remains valid as long as n > 2. Indeed, if we choose pi,..., pn -2 as 
above (all mutually different and ^ p), then if p 7 ^ p, the Ihs reduces to the expression 
± dXP jdp^^ as before, while the rhs vanishes since in that case, p must appear among the 
indices pi,..., pn- 2 ) whereas if p = p, the equation assumes the form 


dXP 

dPi 

which implies 
dXP 


dX’' 


-l/ppl...pn-2 


dp 


dX^ 


u '-pvpi...pn -2 
i 




dX" 


dp\ dp'l 


dXP 

dPi 


(no sum over p and u) , 


(no sum over p, u and k) , 


for mutually different p, u, k. On the other hand, when A^ = l and n = 2, these arguments 
fail, and the only conclusion that can be drawn from equation (I77p is that the following 
divergence must vanish:® 

dXP 

= 0 . (78) 


dp 


Term No. 5: For any choice of indices i,j, p, v and mutually different indices pi,..., Pn-i, 
contracting L^uj^ with the multivector field d^/\ d^ a dp_^ a ... a dp^_^ gives the relation 


d 

i 


dp 


5ZX^ - X° 


dH 

dp'^. 


'-api...pn-l 


_d_ 

dp’'. 


51X^ - 


dH 
F 

i 


dp 


'-crpi...pn-l ’ 


®In the case N = 1, the “internal” index on the position and multimomenturn variables can only assume a 
single fixed value, say i, and so we could in principle just omit it, or else repeat it as often as we like. 
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so that, for any choice of indices i,j, /r, v, p, taking pi,..., pn-i to be 7 ^ p shows that 


dXP dH ^^dX^ dXP dH 
dp'^- ^ dp'' dp'' dpf 

When A^> 1 or n>2, we can use the result of the previous item to conclude that 


(79) 


" ^ dp'< ■ 

Now if to begin with, we fix only the indices z, j and p, we can always choose the other free 
indices u and p in this equation to be equal and 7 ^ p: this reduces the Ihs to the expression 
±dX^ /dp^, while the rhs vanishes. Thus we conclude that 


dX^ 

dpi 


= 0 . 


On the other hand, when N = 1 and n = 2, equation (I79p reduces to a simple statement of 
symmetry:® 

e^’^fsP^ - = 0 . (80) 


dpi dpi dp 


1 dp'; 


Term No. 6 : For any choice of indices i,j, m, v and mutually different indices pi,..., Pn- 2 ) 
contracting with the multivector field d^ a d^ a 9™ a d^^ a ... a dp^_^ gives the 

relation 


dq 




dqi 


'-ucrpi...p„_2 ■ 


As before, this implies 


dXP 

dq' 


= 0 


except when = 1 : in this case, the whole term vanishes identically and no conclusion 
can be drawn. 


To proceed further, we write down the equations obtained from the remaining terms: 
• Term No. 1: 

d /..,dH ..„dH\ d 


9g* V dq^ ^ dp''. 


d 


dxP 


X. , 


(81) 


Term No. 2: 


d 

dpi 


dqi ^ dp'; 


A. 

dx'' 


f) M 


(82) 


Term No. 3: 


A. 

dp; 


Xf + XP ^ 16P X^ - XP ^ 

* dqA dq^ dp'' 

d dXP 2 dx^ 

+ 5 - - 5P 5 ; 


* dx'' 


^ ^ dx^ ' 


(83) 
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• Term No. 4: 


dX^ dX>" 
dq^ dq^ 


(84) 


Assuming that X is projectable to E, we observe that equation (1831) can be integrated directly 
to conclude that 


Xf = H 


dX^ 

dq^ 


Pi 


dX^ 

dq^ 


+ Pi 


dx’^ 


Pi 


dX^ 

dx'^ 


+ T. 


where the are local functions on E which, once again, are independent of the multimomentum 
variables Substituting this relation into equation (1841) . we get 


dY/^ dYj^ 

dq^ dq^ 


which can be solved by setting 


yp 


dY^ 

dq^ 


where the Y!^ are local functions on E which, as before, are independent of the multimomentum 
variables Finally, substituting this expression into equations (f8T]l and (l82]) . we get 


dq^ dx'^ ^ dx^ dx^ dx^ dq^ dpj ^ J 


0 , 


and 

dp^ dx’' ^ dx'^ dx‘' dx‘' dq^ dp'j ^ J 

showing that 


0 , 


H 


dX^ 

dx'^ 


dX^ dYP 


Pi 


dx’' 


+ 


dx’' 


+ 


dH 

dx’^ 


A" 


dH 


dH 


x; = 


Y_ 


where Y_ is a local function on M which is independent of the position variables q^ and multi¬ 
momentum variables Writing Y_ as a divergence. 


dY'^ 

Y = —^ 
dxt^ ’ 

and putting X^ = Y!^ — y!_^, we arrive at equations (I74p and (|75p . 

All that remains to be shown now is the final statement concerning exact hamiltonian vector 
fields. To this end, we apply the exterior derivative to equation (ITHp and subtract the expression 
in equation (f72P : then collecting terms, we get 
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— d{ix ^y) — ix 


+ 

+ 


dxf^ 




i ' o ,, Pi 


dx’- 


-Pi 


dx’^ 


\ dq^ 

dX^ „ ^X^^ \ ^ ^ 


a vfjL \ 

-^H +Xfjd,‘Ad% 


dpj 

r)X^ 

pf dg* A dq^ A d'^x^y 


~*~ dq^ 

BX^ 

+ ^pr dq^AdptAd-x,y. 

Numbering the terms in this equation from 1 to 5, we see that the conditions imposed by the 
fact that X should be exact hamiltonian are the following: 


• Term No. 5: This term vanishes if and only if the coefficients X^ do not depend on the 
variables 

• Term No. 3: Due to the previous condition, this term vanishes if and only the coefficients 
X* do not depend on the variables p^. 

• Term No. 4: This term vanishes if and only if the coefficients X^ do not depend on the 
variables q^, except when X = l: in this case the whole term vanishes identically and no 
conclusion can be drawn. 

• Term No. 2: This term vanishes if and only if the coefficients X^ are defined in terms of 
the coefficients and X* according to equation (174)1 . with X^ = 0. 

• Term No. 1: This term vanishes if and only if equation ()75p is required to hold, with 

X'^ = 0. 


□ 


Remark 2 The classihcation of locally hamiltonian vector fields provided by Theorem [2] is 
not quite complete since it does not cover non-projectable locally hamiltonian vector fields. 
This may not be a reason for great concern since such vector fields are pathological in the 
sense that their flows do not respect any of the bundle structures involved and, perhaps more 
importantly, since such vector fields can only exist in one very special and exceptional case, 
namely when N = 1 and n = 2. Still, it is somewhat annoying that they do not seem to admit 
any reasonable classification. To give an idea of what is involved, consider first the more general 
case of one degree of freedom in any space-time dimension (X = 1, n> 1), where it is common 
practice to omit the “internal” index i on the variables q^ and p7; it is then appropriate to 
redefine the components of X in equation ()7ip . say by writing 


X 


XM 


d 

dx>^ 


~ d 

+ X — + X'^ 
oq 


d 

dpt^ 


(85) 
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Then when n = 2, the arguments presented in the proof above do not allow to conclude that the 
coefficients and X are independent of the multimomentum variables p'^, but only that 


as stated in equation (1751) . 


dpii 


0 , 


dX dX'^ dH 

dpi^ dpi^ dp'^ 


( 86 ) 


(87) 


which, taking into account equation (1781) and using the relation — S'^6^ (with 

any hxed sign convention for ±), is easily shown to be equivalent to equation (|80p . Introducing 
a new function F dehned as 


F 


X - 


dH 

dpf^ ’ 


( 88 ) 


together with the Hessian matrix of the hamiltonian function H, 


H. 


flU 


d‘^H 

dpf^ dp'^ ’ 


(89) 


and its inverse 77^^, it is possible to express all coefficients of X in terms of F and H and their 
partial derivatives up to first order (for T) or second order (for H): 


XM 




dF 

dp'^ 


X 


F - 77'^'^ 


dH dF 

dpt^ dp'^ ’ 


(90) 

(91) 


X^^ 


Hi^u 


dF 


^dx'^ 

_ pjllK, 


+ 


dH dF 
dp^ dq 


d^H 

dq dp'^ 


d‘^H 


d^H 


dx^ dp^ dx^ dp^ 


+ 


dH d‘^H 
dp^ dqdp^ 


dH d^H \ 

dp^ dqdp'^J dp'^ 


(92) 


However, finding the general solution of the entire system seems to be an exceedingly difficult 
task, except if one makes some simplifying assumptions on the hamiltonian 77. One obvious 
choice would be to take 


H = \g^,^{x)p^^p'' + A^(x)p>^ + V{x,q) , 


(93) 


where g represents a Lorentz metric, H is a gauge potential and V is some scalar potential, 
but even in this situation we have not come to a dehnite conclusion. The only case in which 
a complete solution has been found is in the absence of external fields, i.e., when the metric 
tensor g and the scalar potential V are both independent of x whereas the gauge potential A 
vanishes, so M is two-dimensional Minkowski space and g is the standard Minkowski metric ry; 
see [271 Appendix]. 
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6 Covariant phase space 

6.1 Symplectic structure on covariant phase space 

One of the most important properties of the covariant phase space S introduced above (see 
equations (f57|) - (i6n]l ') is that it carries a naturally defined symplectic structure [5l[6l[35] which 
can in fact be derived immediately from the multisymplectic structure on multiphase space |15j . 
Namely, generalizing the prescription of equation ([S} in the sense of using ordinary differential 
forms on multiphase space to produce functional differential forms, rather than just functionals, 
we can define functional canonical 1-forms 2-forms on C, where S is a hypersurface 

in M (typically, when a Lorentz metric is given, a Cauchy surface) and Kj] runs through the 
compact submanifolds of S which are the closure of their interior in S and have smooth boundary 
dKj^, by setting 

JKt. 

for 4> £ Q and S T^Q with X vertical, and 


JKt. 

for (/> e C and £ T^Q with Xi,X 2 vertical. (The same formulas continue to hold if 

we require the vector fields X, Xi, X 2 to be only vertical on the image of 0.) As observed, e.g., 
in [5l[6l[35] and, in the present context, in [15], the restriction of the form to S does not 
depend on the submanifold Ky,, provided that appropriate boundary conditions are imposed:]^ 
this happens because when <() e § and ^Xi'^-: ^X2^ £ T^S, the expression under the integral in 
equation (IMI) is a closed form (called the “symplectic current” in |5|), so that according to Stokes’ 
theorem, its integral over any compact snbmanifold without boundary vanishes. Thus, at least 
formally, covariant phase space becomes a symplectic manifold - albeit an infinite-dimensional 
one; its symplectic form will in what follows be simply denoted by and is explicitly given by 
the formula 

^<I>{^Xi4>,^X24>) = [ 4>*{^X2^x-^<^'h) (96) 

JT. 

where S is any Canchy surface in M and where i^eS and £ T^S, with Xi^X 2 

vertical (or possibly just vertical on the image of (p) and such that supp n supp 5x2^ n s 
is compact. 

6.2 Functional hamiltonian vector fields and Poisson brackets 

The central resnlt obtained in Ref. [T5| can be summarized in the form of two theorems which we 
state explicitly because they form the background for the work reported here. The basic object 
that appears there is the Jacobi operator 3 [</>], obtained by linearizing the De Bonder - Weyl 

®For example, when we compare the integral in equation (I 95 II over two compact submanifolds Jfi.Ei and ^"2,22 
of hypersurfaces Ei and E2 in M whose union U K2,y^2 i® th® boundary dK of a compact submanifold K 

of M, we will get the same result. This also happens when Uif2,E2 i® just part of the boundary of a compact 

submanifold K of M but the remainder, dK \ U K2,t,2)i hu® empty intersection with the intersection of 

the base supports of Sxi 4 ' und 
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operator around a solution cj) € § and whose kernel is precisely the space T^S of solutions of 
the linearized equations of motion, and its causal Green function 

Theorem 3 Given a functional F with temporally compact support on covariant phase space §, 
the functional hamiltonian vector field Xp on S associated to F, as defined by the formula 

fieS, 6(j)e T^S , (97) 

is given by “convolution” of the variational derivative of F (see equation with the causal 

Green function of the corresponding Jacobi operator: 

= f forcf^s. (98) 

Jm o(pJ 

Note that the condition that F should have temporally compact support will guarantee that 
both sides of equation (IU7|) make sense provided we interpret T^S as being the space of solutions 
of the linearized equations of motion of spatially compact support, i.e., we regard it as the 
subspace of the space T^C given by equations (IUT]) - (fMl) where the latter is defined according to 
equation (fT8]l . 

With this statement at hand, it is easy to write down the Poisson bracket of two functionals 
F and G on S: in complete analogy with the formula {f,g} = ix^ixf^ ~ ~dg{Xj) from 
mechanics, it can be defined by 

{F, G}[(/>] = - G'[(f][Xp[(l)]) for (()€§, (99) 

or more explicitly, 

{F,G}[(/>] = - j^d'^x-^[(l)\{x) Xp[4>\{x) for(()GS. (100) 

Combining this expression with that given in Theorem[3l we arrive at the second main conclusion: 


Theorem 4 Given two functionals F and G with temporally compact support on covariant 
phase space §, their Poisson bracket {F, G}, with respect to the symplectic form introduced 
above, is precisely their Peierls - De Witt bracket, given by 

{F,GM = [ d^x [ d^y ^mx) G^^\x,y) -^my) forcfG§. (101) 

Jm Jm <50' 

Note that in view of the regularity conditions imposed to arrive at these results, the previous 
constructions do not apply directly to degenerate systems such as gauge theories: these require 
a separate treatment. 

6.3 The main theorems 

In this subsection, we present the two main theorems of the present paper which, for local 
functionals of the form given by equation ([5]), provide a simple algebraic construction of the 
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functional hamiltonian vector field associated with such a functional and, as a corollary, a simple 
algebraic formula for the Poisson bracket of two such functionals. That such formulas should 
exist is not at all obvions, taking into account that the corresponding formnlas for general 
functionals, as given in Theorem [3] and Theorem 0] above, are essentially analytic: to apply 
them in concrete examples, one needs to solve a system of (linear) partial differential eqnations 
in order to calculate the corresponding causal Green fnnction. Surprisingly, for local fnnctionals 
of the form given by eqnation ([5]), the “convolntion type” integral operator which has this Green 
function as its kernel collapses. 

Theorem 5 Suppose we are given a fiber bundle E (the field configuration bundle) over an 
n-dimensional globally hyperbolic space-time manifold M and a hamiltonian T-L : J®E —> J®E, 
which is a section of extended multiphase space J®E over ordinary multiphase space J®E, 
together with a hamiltonian (n — l)-form f on J®E of spatially compact support such that the 
corresponding hamiltonian vector field X^ on J®E, defined by the formula i^^voy^ = df, is 
projectable to Then given any Cauchy surface T, in M and writing ^ for the local func¬ 

tional on § associated to S and f, as in Definition^ and Xj-^ ^ for the functional hamiltonian 
vector field on S associated to p as in Theorem\^ we have 

= ^Xf4> forfiGS. (102) 

Note that this is a version of Noether’s theorem, i.e., we are dealing with a conservation law, 
for the Ihs of equation (|in2p appears to depend on S and the theorem states that it is equal to 
the rhs, which does not! 

Proof: Rather than analyzing the integral formula (I98|) . we shall show directly that the 

algebraic formula (11021) satisfies all required conditions. First of all, we note that has 

spatially compact snpport because / does. Moroever, it is clear from the results of Section 4 
that equation (|102l) does provide a fnnctional vector field not only on C but also on § since, 
according to equation (l63]l . X^ being locally hamiltonian with respect to ujy and projectable 
to M implies that G T^§> when 0 G S. Therefore, all that needs to be verified is that the 

expression given in equation ()102p satisfies the condition (l97)l . To this end, let Xj denote the 
projection of X^ to M and, for any given (p € S, apply Lemma[I]to construct some projectable 
vector field X'^ on J®E which is ^-related to X^ then the difference Xj- — X^ will be vertical 
on the image of (p, and according to equation ([35]) . equation (I102p becomes 

Inserting this expression into eqnation (I96p . we get, for any dx<P G T^S where X is some vertical 
vector field on J®E of spatially compact support 

^<t>{>^x^j[^],dx4>) = ^<t>{d^^ _^^^fi,dx<p) = j (p*[i^i^^ 

f f o ^) f f 

Now observe that the expression 

^^Recall that according to Theorem [ 2 ] this is automatic if > 1 . 


28 









vanishes due to the assumption that 4> is a solution of the equations of motion. Therefore, 

= J^4'*{ix'^Xf‘^n) ) 

which, according to Proposition [H is equal to 

^ej[^]'^x4> = J^4'*{ixdf) = J^4'*{ix’^Xf^u) ■ 

□ 

An immediate corollary of this theorem is that we can express the Peierls - De Witt bracket 
between two local functionals associated to hamiltonian (n —l)-forms directly in terms of their 
“multisymplectic Poisson bracket”: 

Theorem 6 Suppose we are given a fiber bundle E (the field configuration bundle) over an 
n-dimensional globally hyperbolic space-time manifold M and a hamiltonian T-L : J®E —> J®E, 
which is a section of extended multiphase space J®E over ordinary multiphase space J®E, to¬ 
gether with two hamiltonian (n — l)-forms f and g on J®E of spatially compact support such 
that the corresponding hamiltonian vector fields Xj- and Xg on J®E are projectable to M 
Then given any Cauchy surface S in M, the Peierls - De Witt bracket between the local func¬ 
tionals j and g on § associated to S and f and to S and g, as in Definition{J\ is the local 

functional associated to T, and any of the multisymplectic Poisson brackets {/, g} on J®E that 
can be found in the literature, among them the simple “pseudo-bracket” defined by equation ^ 
as well as the modified bracket defined by equation m- In other words, with any one of these 
choices, we have 

Proof: Combining equation (|99p with equation (I102p from the previous theorem and applying 

equation (131|) from Proposition [T] (which is applicable because the relevant integral need only be 
extended over a compact subset of S such that the base supports of / and g are contained in 
its interior), we obtain, for any 

= “ J^{^*iXfdg) - ixf,M{^*d9)) ■ 

Using that dg = ix we see that the second integral vanishes due to the equations of 
motion dSSD, and we get 

}[</*] = J^4’* (ixg^Xf^H^ = ’ 

where the second equality is obvious if we employ the “pseudo-bracket” of equation ([8]) but 
holds equally well if we employ the modified bracket of equation (fTU|l since, once again, the 
relevant integral need only be extended over a compact subset of S such that the base supports 
of / and g are contained in its interior and then the integral over the additional term can, by 
Stokes’s theorem, be converted to an integral over the boundary of that compact subset and 
hence vanishes. □ 
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7 Conclusions and Outlook 


In this paper we have established a link between the multisymplectic Poisson brackets that have 
been studied by geometers over more than four decades and the covariant functional Poisson 
bracket of classical field theory, commonly known as the Peierls - De Witt bracket. This link is 
based on associating to each differential form on the pertinent multiphase space a certain local 
functional obtained by pulling that form back to space-time via a solution of the equations of 
motion and integrating over some fixed submanifold of space-time of the appropriate dimension, 
considering the result as a functional on covariant phase space (the space of solutions of the 
equations of motion). Here, we have restricted attention to forms of degree n—1, where n is 
the dimension of space-time, which have to be integrated over submanifolds of codimension 1 
(hypersurfaces), but we cannot see any obvious obstruction to extending this kind of analysis to 
forms of other degree. This would be of considerable interest since in physics there appear many 
functionals that are localized on submanifolds of space-time of other dimensions, such as: values 
of observable fields at space-time points (dimension 0), Wilson loops (traces of parallel transport 
operators around loops) in gauge theories (dimension 1), electromagnetic field strength tensors 
and curvature tensors (dimension 2), etc. 

The overall picture that emerges is that the correct approach to the concept of observables 
in classical held theory is to regard them as smooth functionals on covariant phase space. 
As is well known, covariant phase space is, at least formally and for nondegenerate systems, 
an (inhnite-dimensional) symplectic manifold, so the space of all such functionals constitutes 
a Poisson algebra, and that is what we are referring to when we speak about the “algebra of 
observables” in classical held theory. Of course, this algebra is huge, and the construction of local 
functionals on covariant phase space from differential forms on multiphase space, as employed in 
this paper, is merely a device for producing special (and quite small) classes of such observables. 
But the reduction of the algebraic structure at the level of such functionals to some corre¬ 
sponding algebraic structure for the generating differential forms is highly problematic: in fact, 
there is no reason to expect that there might exist any product or bracket between differential 
forms on multiphase space capable of reproducing the standard product or bracket between the 
corresponding functionals on covariant phase space. One possible obstacle is that any such pre¬ 
scription would most likely be highly ambiguous since a crucial piece of information is missing: 
after all, the functional does not only depend on the differential form which (after pull-back) is 
being integrated but also on the submanifold over which one integrates! One way out would be 
to restrict to functionals defined by some fixed submanifold and hope that the resulting algebraic 
structure does not depend on which submanifold (within a certain given class) is chosen: that 
is what we have done in this paper when reducing the covariant Peierls - De Witt bracket to a 
multisymplectic bracket. But the fact that this actually works is to a certain extent a miracle 
which cannot be expected to happen in general, since we will very likely be forced into admit¬ 
ting functionals defined by integration over different submanifolds, including submanifolds of 
different dimension. For example, this happens as soon as we want to include differential forms 
of different degrees and/or explore the existence of a relation between the product of functionals 
and the exterior product of forms, or some modified form thereof. In particular, these argu¬ 
ments show why the notorious absence of a decent associative product in the multisymplectic 
formalism should come as no surprise: it merely expresses the fact that functionals of the form 
given by equation ([5]) do not form a subalgebra. 
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At any rate, it is a highly interesting question what kind of algebraic structure on what kind 
of spaces of differential forms (or pairs of submanifolds and differential forms) will ultimately 
result from the functional approach advocated in this paper. These and similar questions are 
presently under investigation and will be reported in a future publication. 
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